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ABSTRACT 

Compact radio emission provides a reliable method for the detection of low luminosity young stellar objects 
\ (YSOs), and is particularly useful for detecting the earliest stages of protostellar evolution where the source 

psj . itself may still be heavily embedded in its natal dust envelope. For such Class and Class I objects the dominant 

radio emission mechanism is expected to be free-free, however unlike massive YSOs the way in which this 
O I radio emission is produced remains a subject of debate. As larger samples of radio YSOs become available the 

relationship between the radio luminosity of the Class O/I population and their wider global properties is now 
being clarified. Furthermore, the broader scientific applications of such samples are also becoming increasingly 
apparent. These improved constraints on the nature of the radio emission from YSOs are now contributing to 
our understanding of not only the evolutionary physics of protostars themselves but also their wider impact on 
their surroundings. Here we discuss the physics of the radio emission, the emerging relationship between this 
' emission and other properties of YSOs and some of the applications for studies exploiting this emission. 

. Subject headings: radiation mechanisms:general - stars: formation - ISM:clouds 

^' 

1- INTRODUCTION 

Q ■ The long wavelength radio emission from low-mass young stellar objects (YSOs; Class & Class I) is capable of escaping 
] from the high column density dust envelopes which surround these deeply embedded young protostars, making such objects 
c/3 . detectable in the radio band even when opacity effects conceal them at infrared frequencies. Although the thermal dust spectrum 
" falls off steeply in intensity at longer wavelengths, young and low-luminosity protostars are known to produce additional radio 
emission. This radio emission has been observed to possess a rising spectrum, indicating that it occurs as a consequence of 
I . free-free radiation from an ionized plasma. Although the physics behind the ionization around high-mass stars is reasonably 
<^ • well understood as arising from photoionization due to the strong UV flux from such objects, in the low-mass case a number 
Qv^ of mechanisms have been proposed to produce the ionization responsible for the observed free-free emission. Distinguishing 
ly-j . between these possible mechanisms is not trivial from currently available data. Where a high enough ionizing flux is present, 
■ generally in later type T Tauri stars, photoionization can also support an embedded Hll region (Churchwell 1990); a fully ionized 
[ — ■ stellar wind, associated more often with more evolved stars than very young YSOs, would produce an observable radio signal 
(Panagia & Felli 1975; Curiel et al. 1989), as could a partially or fully ionized, collimated outflow (Reynolds 1986); the accretion 
shock on the surface of protostellar discs may heat and ionize infalling gas (Winkler & Newman 1980; Cassen & Moosman 
1981); molecular outflows from young, low-mass protostars and the neutral winds thought to initiate them can generate free- 
free emission through shock ionization as they impact on the surrounding envelope (Curiel, Canto & Rodriguez 1987; 1989; 
Rodriguez & Reipurth 1996), a mechanism which was refined by Ghavamian & Hartigan (1998). 

It is often the last of these mechanisms which is considered to be most plausible and indeed, where measured, the outflow force 
of many protostellar jets has been found to be energetically viable to explain the observed cm-wave radio emission (Anglada 
^ , 1995). This theory is also supported by very high angular resolution observations of the radio emission from protostars (Simon 
' et al. 1983; Anglada 1996), which often find that the radio emission is elongated along the direction of the molecular outflows. 
■ - ' However, as increasingly sensitive molecular observations reveal weaker and weaker outflows (e.g. Bourke et al. 2005; Crapsi 
et al. 2005) the measured momentum flux is also found to be insufficient to explain the total radio emission in a rising number 
of objects (e.g. Scaife et al. 201 la) and it is possible that a combination of the above mechanisms are infact responsible for the 
observed radio emission. 

As mentioned above, it is generally assumed that the dominant mechanism giving rise to continuum radio emission from YSOs 
is thermal bremsstrahlung, or free-free, emission. This mechanism operates when charged particles, typically electrons, are 
accelerated by encountering another charged particle. Consequently radio free-free emission is expected to manifest from any 
environment populated by an ionized plasma. The distribution and energy cut-off for scattering in such a plasma is characterized 
by the Gaunt factor, (gff), which allows one to work out the (ratio of) absorption coefficient, Ky, for the plasma.Integration of 
this quantity along the line of sight, / KyM, gives the optical path length (depth) for the free-free emission, Ty. Optical depth 
is a strong function of electron density, generally parameterized through the emission measure of the plasma, EM = J NeNiM, 
and a weaker function of electron temperature. Once the optical depth of the emitting plasma has been determined the brightness 
temperature, Tb, of the resulting emission is given by the product, 

rb = 7;(i-e-^^). (1) 

In the Rayleigh- Jeans region of the spectrum the resulting flux density is then given by Sv^2kjiTb/X^Q.. 
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Where the plasma giving rise to the free-free emission is reasonably uniform in density this leads to a characteristic radio 
spectrum which has two components delineated by the frequency at which the optical depth equals unity (Ty = 1), marking 
the transition from optically thick behaviour to optically thin. In the optically thick regime the optical depth term becomes 
approximately unity and the flux density spectrum rises as V^; in the optically thin regime the optical depth term becomes « Ty 
and the frequency dependence of Ty cancels with the exphcit frequency dependence in the Raylegh-Jeans approximation, resulting 
in a flux density spectrum which varies as v^'* '. 

Key to free-free emission is the presence of a hot, ionized plasma. For massive YSOs photoionization of surrounding neutral 
gas by the ionizing UV flux from the central object is high enough to account for the observed radio emission. For low mass, 
low luminosity objects this is not the case. Fig. [1] shows the expected radio luminosity from the Lyman continuum flux of a 
given bolometric luminosity 'zero-age-main-sequence' (ZAMS) star This comparison was first shown by Anglada (1995) and 
following that work we adopt Lyman continuum fluxes from Thompson (1984) for a given range of bolometric luminosities. The 
radio luminosity at a given frequency can then be derived using 

(4^) -2.08X 10--f^) i-^y-' fiV", (2) 
VmJykpcV Vctss-iy VGHz/ \K J ' 

which has been adapted from Wilcots (1994). Here we assume ~ 10^ K and update the measured data to include objects at 
lower luminosities from the AMI-LA radio continuum follow-up (AMI Consortium: Scaife et al. 201 1 a;b; 2012) to the Spitzer 
c2d survey (Evans et al. 2005). 

It is immediately obvious from Fig.[T]that the radio luminosity of low bolometric luminosity YSOs is significantly in excess of 
that expected from photoionization, and that this trend continues to very low (Lboi ^ 1 Lq) luminosities where the discrepancy 
reaches > 10 orders of magnitude. 

The alternative explanations for the radio emission associated with these low luminosity YSOs as listed above generally invoke 
the action of powerful molecular outflows and stellar winds associated with such objects. The presence of these energetic 
processes provides an alternate source for the ionization required to explain the observed free-free emission. 

In these cases the ionized plasma generally cannot be well approximated by uniform optical depth but instead has regions of 
varying density, resulting in a partially opaque plasma. The spectra from such regions have behaviour that depends not only on 
the thermodynamic conditions, but also on the geometry of the ionized gas. 

The canonical geometry for a stellar wind is that of a spherical region of ionized, isothermal' plasma with a radial density 
gradient, ne(r) = no('"/'"o)^ (Panagia & Felli 1975), where «o is the electron density at the inner boundary ro, which marks the 
stellar radius. For such an emission region the optical depth through the plasma varies as a function of projected radius on the 
sky. 

For a generalised outflow where the physical conditions are parameterized by different radial dependencies such that Te{r) = 
(r/ro)'^^ , T(r) = (r/ro)^^ etc the spectral index of the radio emission would have the form 

2 1 

a = 2 + — (1 + e + ^r) (3) 

with q-i = £ + 2qx + 2q„ — 1.35qT (Reynolds 1986), where T denotes optical depth, x denotes ionization, n denotes density, T^. 
denotes temperature and £ relates to the opening angle of the jet such that w — wo{r/roY with w representing the radius of the 
jet perpendicular to the direction of outflow. For a range of physical situations this index varies with values —0.1 < a < 1.1, as 
listed in Table 1 of Reynolds (1987). 

Opening angles for outflows from young stellar objects (YSOs) vary with evolution. For example the earliest stage of YSO 
(Class 0) have highly collimated outflows which become less collimated as they evolve through Class I and, where an outflow 
remains, leave a wide angle outflow from Class II objects. Typically, however, one may assume that Oq < 0.5 rad. Outflows 
which are collimated (e < 1) will have shallower radio spectral indices, typically a = 0.25, although pressure confined outflows 
can exhibit much steeper indices approaching unity. Such steep indices can also be produced by recombining or accelerating 
outflows where the indices may even exceed unity. A conical, fully ionized, isothermal wind will have the same spectral index as 
the canonical spherical stellar wind of Panagia & FelH (1975): a ~ 0.6. This is assumed to be the typical spectral index for radio 
emission from YSOs, and indeed the radio spectral indices from samples of YSOs, where measurable, typically have a mean 
value of 0.6 (e.g. Scaife et al. 2012). 

2. CHARACTERISTIC CORRELATIONS 

In most cases the number of free parameters in detailed theoretical models of partially optically thin free-free emission far ex- 
ceed the physical constraints available from observational data. Consequently it is typical to express the expected radio luminosity 
as a function of alternate observables. 

It can be seen from Fig. [T|that there is a strong correlation between radio and bolometric luminosity (r = 0.77; Scaife et al. 
201 la), characterized as 

logL,ad[mJy kpc2] ^ -(1.74±0.18) + (0.51 ±O.26)logLboi[L0], (4) 
at a frequency of 16 GHz (Scaife et al. 2011b). An even stronger correlation (r = 0.89; Scaife et al. 2011b) is evident when 

' The polytropic case for the same geometry was considered by Chiuderi & Torricelli Ciamponi (1978), whose derivation pre-empted the analysis of Reynolds 
(1986) but in a less general form. 
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Fig. 1 . — Observed radio luminosity as a function of bolometric luminosity, (a) The solid line shows the radio emission expected from Lyman continuum 
radiation from a ZAMS stai; see text for details. Data points are from the AMI-LA (filled circles; AMI Consortiuin: Scaife et al. 201 la; b; 2012) and from the 
compilation of Anglada (1995; crosses), (b) L^j — Lbol coirelation showing compiled data froin the literature norinalized to 5 GHz assuming a global spectral 
index of a = —0.6. 



considering the relationship between radio luminosity and envelope mass, 

logLrad[mJy kpc^] = -(2.23 ±0.65) + (0.68 ±0.62) logMenv[Mc: 



(5) 



(Scaife et al. 201 lb). Similarly to the correlations of outflow force, Fqo, with bolometric luminosity and envelope mass charac- 
terized by Bontemps et al. (1996) there is a stronger dependency on evolutionary state indicated in the relationship between radio 
luminosity and envelope mass, compared to that with bolometric luminosity. The power-law index of the relationship Z^ad °= ^env 
is dominated by the presence of Class I objects, which deviate from the general trend (Scaife et al. 201 lb). When considering 
Class objects alone the power-law index is significantly increased and in both cases is in fact consistent with that expected from 
models of both spherical (Panagia & Felli 1975) and collimated (Reynolds 1986) stellar winds, which predict that 
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as for simple models of protostellar evolution the envelope mass is expected to be proportional to the mass loss rate (Bontemps 
et al. 1996; Fuller & Ladd 2002; Arce & Sargent 2006). 

Anglada (1995) suggested that there should also be a correlation of radio luminosity with outflow force, as might be expected 
from the shock ionization models of Curiel et al (1989). Where molecular outflows from young, low-mass protostars and the 
neutral winds thought to initiate them generate free-free emission through shock ionization as they impact on the surrounding 
envelope (Curiel, Canto & Rodriguez 1987; 1989; Rodrguez & Reipurth 1996) the free-free flux density would be proportional to 
the outflow force (also called momentum flux), Fco, which is equivalent to the rate of outflow momentum and is often calculated 

asFcO = PCo/fdyn- 

At a wavelength of 2 cm the expected flux density from such shock ionization can be computed as 
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(Scaife et al. 2010; Anglada 1996), where rj — {D./4n) is the fraction of the stellar wind being shocked and the factor £, (t) = 
t/(1 — e^^) allows for the fact that the radio emission may not be optically thin (Anglada et al. 1998), although this dependence 
on the optical depth, T, is very weak. 
Fig. |2|a) shows the distribution of observed radio luminosity with outflow force for data compiled from the literature. A 

weak positive correlation is present in these data with r = 0.62 and a best fitting regression of Fqo ^ ^j-a(j ' consistent with 
that determined by Anglada (1995) for the original dataset. Removing the luminosity dependence from these data by using the 
relationship determined by Bontemps et al (1996), Fqo °= Lboi, and that shown in Eq.|4l appears by eye to improve the correlation, 
see Fig.|2|b) but in fact results in an even poorer Pearson co-efficient (r = 0.36). In both cases the regression is consistent with 
that expected for the shock ionization model, but the correlation of radio luminosity with outflow force remains very weak with 
substantial scatter 

There are a number of effects, resulting from the different methods used to calculate outflow force, that will add scatter to these 
distributions. While some estimates are derived from the total momentum observable in the outflow lobes, divided by a measure 
of the outflow dynamical time, others use a measure of the force derived within an annulus around the outflow origin (Bontemps 
et al. 1996). This later provides a better estimate of the force for outflows which have been only partially mapped. Further 
variations arise from whether the maximum velocity in an outflow is used to compute the dynamical time or some average/typical 
velocity. Flux loss in interferometric studies, high optical depths of CO emission in the outflow line wings and inclination effects 
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Fig. 2. — Observed radio luminosity as a function of outflow force, (a) Compiled data from the literature are shown normalized to 8.35 GHz assuming a global 
spectral index of a = —0.6. The minimum outflow force required to produce a given radio luminosity is indicated for plasma with electron temperatures of 
3000 K (dashed hne) and lO'* K (solid line), (b) Fcoc/ibol (dimensionless) versus M^^^/lJ^^-^ for all objects in the listed data sets that have available values for all 
three physical parameters. 

will all also reduce the observed outflow velocities and tend to reduce the momentum actually observed. Many studies follow 
the lead of Bontemps et al. (1996) and increase their outflow forces by a factor of 10 to account for optical depth and inclination 
effects (e.g. Visser et al. 2002). Optical depth effects can also be corrected for using observations of the rarer isotopologues 
of CO (e.g. L1221-IRS3; Umemoto et al. 1991), and some studies note that the inclination corrections will be minimal (e.g. 
L1014-IRS; Bourke et al. 2005). 

In addition to the lack of a strong correlation between radio luminosity and outflow force, a further concern from these data 
is the number of objects which appear to have outflow force values below the minimum level required to explain their radio 
luminosity via the shock ionization model. This level is indicated by a solid line in Fig.|2a). Such sources are evident even in the 
initial compilation of Anglada (1995) and the population has become more pronounced as an increasing number of objects with 
weaker outflows have been detected (Scaife et al. 201 lb; 2012). This discrepancy is important as it may indicate that the shock 
ionization model for radio emission from YSOs is not the dominant mechanism operating in these objects. 

However, the observed discrepancy may also be an effect of the method used for calculating outflow force: it was noted by 
Curtis et al. (2010) that the dynamical time may not be suitable as a divisor for calculating the outflow force. Dynamical time, 
defined as Td = Liobe/t^'max, is the time required for the bow shock of an outflow to travel the distance of that outflow at the 
maximum velocity observed in the outflow. This quantity is not a good measure of the age of the outflow. It has been suggested 
both to underestimate true outflow ages in evolved sources (Parker, Padman & Scott 1991) and to overestimate outflow ages in 
young Class sources (Downes & Cabrit 2007). A better tracer of outflow age has been suggested to come from the bolometric 
temperature of a source (Curtis et al. 2010), which is expected to increase as a protostar consumes its surrounding dust envelope 
and advances towards the main sequence. 

The age of the driving source, and hence outflow age, can be estimated from the bolometric temperature using the relation. 



logfYSoN - (-0.9±0.6)+2.41og7boi[K] 



(8) 



(Ladd, Fuller & Deane 1998). Examining a subsample of the data presented in Fig.|2|a), which has information available to 
calculate both the typically used outflow force using the dynamical time , Fqo = pco/'^d, and the outflow force calculated using 
the outflow age, F^q — pco/tyso, it can be seen that objects which are inconsistent with the Curiel et al. (1989) shock ionization 
model using the first method, become consistent when the second method is applied. This distinction is illustrated in Fig. [3] 
The one object which remains inconsistent with the shock ionization model is the only Class I source in this sample (circled), 
suggesting that for Class I objects an additional radio emission mechanism may be present. In addition, although only one Class I 
object is present in this sample we note that the object currently known to be most discrepant with the shock ionization model in 
this regard (L1014-IRS1; Shirley et al. 2007; Scaife et al. 201 la), shown as the most low-lying point in Fig.|2a) but not present 
in the subsample, is also Class I. 

3. RADIO SELECTION FUNCTIONS 

In spite of the strong correlation between radio and bolometric luminosity, this relationship is defined by protostellar objects 
which are detected in the radio. The distribution of non-detections of protostellar objects does not seem to be a function of 
the radio luminosity detection limit when considered against bolometric luminosity, but is instead better traced by the radio 
luminosity detection limit considered against envelope mass (Scaife et al. 2012). Indeed those objects which are not detected 
in the radio are sources that have a high ratio of bolometric luminosity to envelope mass. Such a ratio is commonly used as an 
evolutionary tracer amongst YSOs; in Fig.|4]we plot bolometric luminosity against envelope mass for a sample of objects from 
the Perseus molecular cloud (Hatchell et al. 2007), the delineating nature of Menv/iboi is shown as a red dashed line and the 
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Fig . 3 . — Data for a subsample of objects is shown with outflow forces calculated in different ways. Fqo = pco / uses the traditional method with dynamical 
time as the divisor, whereas Fqq = pco /'yso uses the alternate method with the age of the driving source as the divisor, see text for details. All objects are 
Class 0, with the exception of the ringed source which is Class I. The minimum outflow force required to explain the observed radio emission using the shock 
ionization model of Curiel et al. (1989) is indicated by a dashed line. 

ratio determined by Bontemps et al. (1996) of Menv/^bo^ ^ ^^^^ dashed line. An evolutionary distinction in radio detection 
rates with a ratio of ^ Mgnv/iboi is consistent with the overall radio detection rates for Class and Class I objects, which are 75 
and 25 per cent, respectively (Scaife et al. 2012). These detection ratios are common across multiple different star formation 
regions and, considered in light of the mass-luminosity ratio trend described previously, strongly suggest that radio detections 
are biased towards younger embedded objects. A consequence of this is the implication that blind radio surveys of star forming 
regions will preferentially detect those objects which have high ratios of envelope mass to bolometric luminosity. This makes 
radio surveys of such regions highly complementary to surveys in the (mid-)infared, as those data will be biased against such 
objects due to opacity effects. However, the very low detection rate for Class I objects, combined with the fact that Class I objects 
exhibit a degree of variability in their radio emission which is likely to arise from alternative emission mechanisms, makes this 
evolutionary behaviour difficult to quantitatively demonstrate at present due to small sample sizes. 

4. RADIO EMISSION AND THE EVOLUTIONARY SEQUENCE 

The population of protostellar objects discussed here encompasses the earliest stages of their evolution: Class and Class I. 
More evolved pre-main sequence (PMS) stars. Class II/III, are also known to exhibit radio emission (see e.g. Feigelson & 
Montmerle 1999). For these later objects the radio emission is dominated by non-thermal processes associated with magneto- 
spheric/coronal activity and the radio emission is known to be strongly correlated with X-ray emission (Giidel & Benz 1993). 
This correlation is thought to be due to a common underlying energy reservoir being responsible for both plasma heating that 
produces thermal X-ray emission and particle acceleration in the presence of magnetic fields, which prodices non-thermal radio 
emission. Establishing whether such a radio/X-ray relation holds for younger protostellar objects is complicated by a number of 
factors. Two factors which cause immediate issues are, firstly, that the thermal emission, which is observed to dominate the radio 
spectrum for Class O/I objects, quickly becomes opaque to non-thermal emission (Andre 1987) and secondly, that the high degree 
of variability in the emission from Class I/II objects requires simultaneous X-ray and radio observations to be taken. Even taking 
these factors into account, a recent study of the radio/X-ray correlation for younger protostellar objects showed that the recov- 
ered data are inconsistent with the extrapolated Giidel-Benz relation by more than an order of magnitude (Forbrich et al. 201 1). 
However, these results are not yet definitive as such studies are limited to small samples of objects, both due to observational 
constraints as well high non-detection rates for as yet unknown astrophysical reasons, as well as an uncertainty in the underlying 
mechanism of the detected radio emission. 

The differing radio detection rates for Class and Class I objects noted in §|3]may indicate a transition between underlying 
emission mechanisms during the Class I stage, consistent with the expectation that some fraction of thermal emission may also 
be expected into the Class II stage. This remnant of thermal activity in the case of Class II objects is perhaps associated with 
ionization due to a surface shock (e.g. Cassen & Moosman 1981) or ionized plasma in a disk wind (Pascucci et al. 2012). A 
combination of mechanisms in the case of Class I objects is also indicated by the polarization signature seen in case of L1014 
(Shirley et al. 2007) during flare activity, which can only have arisen from a non-thermal mechanism. The radio detection rates 
for Class O/I objects are inverted in the X-ray, for example in the Serpens and NGC 1333 regions no Class objects were detected 
using Chandra, whilst 23% of Class I objects were detected (Winston et al. 2010). It is unclear at this stage whether this effect 
is due entirely to extinction/absorption of soft X-rays or if it represents a fundamental underlying trend. Future work in this 
area, examining the relationship between radio and X-ray emission, will be instrumental in distinguishing the underlying radio 
mechanisms and their distribution across protostellar class. 

5. RADIO EMISSION AS A TOOL 

Historically the study of radio emission from low mass young stellar objects (YSOs; Class & Class I) has concentrated 
on characterizing the radio emission itself and the underlying emission mechanism. However, recent advances in a number of 
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Fig. 4. — Relationship between envelope mass and bolometric luminosity. Data are shown for the Perseus molecular cloud from Hatchell et al. (2007) with 
Class objects shown in blue and Class I objects shown in green. The envelope mass detection threshold for these data is indicated by a light blue shaded area. 

Other areas have shown that understanding, quantifying and predicting the radio emission from protostellar objects is increasingly 
important for a variety of reasons, including, but not limited to, the investigation and characterization of the anomalous microwave 
emission attributed to spinning dust, tackling the luminosity problem indicative of an incomplete understanding of protostellar 
accretion history, and exploring the long wavelength tail of the greybody distribution in order to understand dust grain evolution 
and planet formation. 

The first of these, the observational study of the anomalous microwave emission attributed to spinning dust, has largely con- 
centrated on arcminute scale emission since it is primarily targeted by CMB telescopes, which consider it a source of foreground 
contamination. Although dark and molecular clouds have been identified as the source of such anomalous emission (e.g. L1622, 
Cassasus et al. 2006; L675 & LI 111, Scaife et al. 2010) these large-scale observations have largely neglected the known 
small-scale radio emission detected from a significant fraction of protostars (Rodriguez et al. 1989; Anglada 1995). As such 
investigations become more detailed it is necessary to quantify correctly this emission. Firstly, in order to investigate any connec- 
tion between the presence of active star formation within a dark cloud and the identification of a spinning dust component. The 
observational uncertainties surrounding the identification of heavily embedded, low luminosity YSOs have lead to the proposition 
that many clouds currently considered to be starless, including a number not previously thought to even be undergoing collapse, 
do in fact host star formation. This possibility creates complications when investigating the correlation of those clouds hosting 
star formation with those which are associated with anomalous microwave emission, as highlighted by the case of L675 (Scaife 
et al. 2010) and this remains an important open question. Secondly, in marginal cases, identifying radio emission from protostars 
internal to the cloud is important in order to avoid confusing their emission with any excess that may be attributed to spinning 
dust. 

Most radio protostar searches and surveys (e.g. Andre, Montmerle & Feigelson 1987; Anglada 1995; Stamatellos et al. 2007) 
are conducted at 3.6 and 6 cm wavelengths; however, the spectrum of emission due to spinning dust is thought to peak at shorter 
wavelengths of 1 — 2 cm (Draine & Lazarian 1998) and this is the wavelength range where most spinning dust observations are 
made. Although historically the number of observations made at 2 cm (Rodriguez & Canto 1983; Bieging, Cohen & Schwartz 
1984; Pravdo et al. 1985; Rodriguez & Reipurth 1998; Andre, Motte & Bacmann 1999) were limited, particularly in the case of 
low-luminosity protostars and Class objects, a substantial increase in the number of such objects observed has recently been 
made through the Spitzer follow-up surveys of the AMI telescope at 16 GHz (Scaife et al. 2011 a;b;2012). These surveys were 
initiated to quantify the radio contribution of, and correlation with, YSOs to the microwave emission of dark clouds following 
the detection of a new radio YSO in L675 (Scaife et al. 2010). 

The strong correlation between radio and bolometric luminosity first noted by Anglada (1995) has been shown by recent 
surveys to hold towards lower luminosities (Scaife et al. 2012), see Fig. [T| Although historically radio observations of protostars 
have been skewed towards higher luminosity objects due to observational constraints, improvements in sensitivity have allowed 
this distribution to tend more towards that seen at infrared wavelengths and predicted by theoretical models of YSO evolution (e.g. 
Dunham & Vorobyov 2012). One consequence of this result is that predicting the contribution of protostars, where identified, 
to the overall microwave emission of star-forming dark clouds may be done in principle from their bolometric luminosity in the 
infrared. A further area where such a relation may prove important is in determining the abundance of low luminosity protostars. 

Measurements of the luminosity distribution of protostars from the Spitzer Space Telescope (Dunham et al. 2008; Evans 
et al. 2009) have aggravated the 'luminosity problem' (Kenyon et al. 1990). This problem exhibits as an increasing number 
of protostars at lower luminosities down to an internal luminosity of Lint — 0.1 Lp, below which the numbers start to decline. 
However, simple star formation models, such as the spherical accretion model of Shu (1977), predict that a low-mass source on 
the stellar/brown dwarf boundary {M ~ 0.08 Mp) should have an internal luminosity Lint — 1.6Lo from accretion alone (Evans 
et al. 2009). Non-steady accretion, starting in the earliest protostellar stages, is currently the best solution to this discrepancy 
(Kenyon & Hartmann 1995; Young & Evans 2005; Enoch et al. 2007). The luminosity problem is most difficult to rectify in very 
low luminosity objects (VeLLOs; Young et al. 2004; Dunham et al. 2008) with extreme luminosities Lint < 0.1 Lq. 
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Recent predictions of the luminosity distribution expected from non-steady accretion scenarios (e.g. Dunham & Vorobyov 
2012) have instead predicted an over-abundance of low luminosity objects: the so-called 'reverse luminosity problem'. These 
results emphasize the point that understanding the true accretion mechanism depends on observationally constraining the shape 
of the luminosity distribution at low luminosities. The existing sample of VeLLO sources, although expanded by Spitzer, is in no 
way complete. VeLLOs are difficult to confirm as protostars in the infra-red due to their low luminosity and embedded nature. It 
is also the case that VeLLOs are often found in cores which are not only assumed to be starless, but which were also not believed 
to be approaching collapse (Bourke et al. 2006). Nevertheless, identifying a complete sample of these low luminosity embedded 
protostars is vital for understanding low mass star formation. Radio emission from heavily embedded objects provides a useful 
tool for establishing the presence of protostellar activity in such cases as it is not subject to the same degree of opacity effect as 
the shorter wavelength infra-red emission. 

6. CONCLUSIONS 

The long wavelength radio emission from YSOs provides a novel tool for confirming the population of low luminosity embed- 
ded protostars. Increased samples of objects have extended the known correlation of radio luminosity with bolometric luminosity 
by over an order of magnitude since it was first noted (Anglada 1995), and this correlation now includes a substantial number 
of the newly discovered VeLLO objects. A further correlation with envelope mass has also been identified. Potential problems 
from sources with small outflow force values for the generally accepted shock ionization model of free-free emission have been 
shown to disappear when outflow force is calculated using the age of the driving sources as opposed to the dynamical time. In 
the case of Class I objects this discrepancy is more difficult to rectify and, when coupled with a lower detection rate, this suggests 
that although shock ionization may be the dominant emission source for Class objects an alternate mechanism may contribute 
substantially to the emission of Class I sources. 

These improved constraints on the nature of the radio emission from YSOs are now contributing to our understanding of 
not only the evolutionary physics of protostars themselves but also their wider impact on their surroundings through e.g. their 
connection to the anomalous microwave emission attributed to spinning dust. As further improved observational samples become 
available these constraints can only improve and it is doubtless that even more applications will be revealed. 
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